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ABSTRACT 

A  computer  simulation  of  an  aircraft  making  a  low  level  penetra- 
tion over  hostile  terrain  is  described.     The  model  discussed  is  intended 
as  a  tool  for  the  determination  of  appropriate  defensive  tactics  for  air- 
craft upon  receiving  fire  from  conventional  ground  weapons.    In  any  com- 
puter model  of  a  ground  to  air  combat  environment  it  is  essential  to  have 
a  realistic  simulation  of  the  terrain.    The  method  used  in  this  simulation 
is  discussed  and  shown  to  be  extremely  accurate. 
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1 .         Introduction 

PENAIR  is  a  time  step  computer  simulation  of  an  aircraft  flight  over 
hostile  terrain.     It  is  an  analytic  computer  war  game  in  which  one  aircraft 
and  a  maximum  of  two  hundred  weapons  are  the  principles.     Any  piece  of 
actual  terrain  may  be  simulated  using  a  contour  map  of  the    area  to  obtain 

the  input  parameters  . 

The  simulation  is  intended   as  a  tool  for  the  determination  of  appropri- 
ate defensive  tactics  for  an  aircraft  upon  receiving  fire  from  conventional 
ground  weapons.    Using  the  characteristics  of  the  aircraft  and  of  the  wea- 
pons the  model  determines  the  survival  probability  of  the  aircraft  on  any 
of  several  flightpaths.    There  are  no  restrictions  on  the  aircraft's  speed, 
altitude  or  other  characteristics  and    therefore  this  model  is  equally  suita- 
ble for  helicopters  or  jets. 

The  aircraft  flightpath  can  be  predetermined  or  it  can  be  generated 
by  the  model.    Upon  receiving  fire  the  aircraft  can  be  programmed  to  exe- 
cute any  of  18  evasive  maneuvers  .    A  comparison  of  the  survival  probabili- 
ties can  then  be  made  thus  giving  an  indication  of  the  appropriate  defensive 
maneuver  to  employ. 

The  terrain  simulation  and  generated  flightpaths  appear  to  be  ex- 
tremely realistic.    The  terrain  is  simulated  by  fitting,  in  the  least  squares 
sense,  a  polynomial  function  of  two  independent  variables  (x,y)  to  values 
of  a  dependent  variable  specified  at  points  on  a  rectangular  grid  in  the 
plane  of  the  independent  variables .    This  is  accomplished  by  a  separate 
computer  program,  TERRAIN,  which  generates  the  necessary  inputs  for 
PENAIR. 

2  .       Design  of  the  Model 

PENAIR  is  a  computer  simulation  of  a  single  aircraft  making  a  low 
level  penetration  over  terrain  containing  hostile  small  arms  and  con- 
ventional anti-aircraft  weapons  .    Any  actual  or  fictitious  terrain  may  be 


simulated.    No  offensive  aircraft  weapons  are  considered  although  18  dif- 
ferent defensive  maneuvers  such  as  increasing  speed  and  diving  for  the 
deck  may  be  performed.    Ten  types  of  ground  weapons  may  be  used  with 
a  maximum  of  20  weapons  per  type.    The  program  computes  the  aircraft's 
time  in  view  of  each  weapon,  the  probability  of  the  aircraft  surviving  each 
weapon  and  the  aircraft's  total  survival  probability.    The  simulation  runs 
in  discrete  time  steps  using  any  time  interval  desired.     The  program  does 
not  utilize  Monte  Carlo  methods;  aircraft  survival  probability  for  a  speci- 
fied flightpath  over  a  specified  terrain  is  computed  in  a  single  run.     Flow 
chart  1.0  of  Appendix  Ir  page  30  ,  indicates  the  general  logic  of  the  PENAIR 
program . 

PENAIR  is  written  in  FORTRAN-63  for  the  CDC  1604  computer.    The 
PENAIR  program  is  included  as  Appendix  II.    A  separate  program,  TERRAIN, 
is  utilized  to  supply  inputs  to  PENAIR  for  terrain  simulation.     The  terrain 
input  program  supplies  the  necessary  data  on  punched  cards  for  the  PENAIR 
program.    The  TERRAIN  program  is  also  programmed  in  FORTRAN-63  for  the 
CDC  1604  and  is  included  as  Appendix  rv. 

The  terrain  used  in  the  PENAIR  program  may  be  divided  into  as  many 
as  12  separate  squares  for  simulation.    Each  of  the  12  individual  squares 
is  approximated  by  a  separate  polynomial.    The  12  squares  may  be  of  any 
common  size.    Figure  1  illustrates  the  arrangement  of  these  squares.  The 
coordinate  system  has  been  established  such  that  the  lower  left  corner  is 
the  origin  with  increasing  x  to  the  right  and  increasing  y  upward.     It  is 
possible  to  use  actual  military  grid  coordinates  since  the  origin  may 
take  any  desired  x  and  y  values . 

A  special  PENAIR  output  is  a  set  of  punched  cards  to  be  used  as 
input  to  a  graph  plot  routine  for  an  off-line  plotter.    This  routine  then  draws 
the  actual  flightpath  over  the  simulated  terrain.    Output  from  this  routine 
is  included  in  Section  6. 
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FIGURE  1 

THE  ARRANGEMENT  OF  THE  INDIVIDUAL 
SQUARES  OF  TERRAIN  SIMULATED  IN 
THE  PENAIR  PROGRAM 


3 .    Aircraft  Movement 

The  PENAIR  model  has  provisions  for  three  types  of  simulated  aircraft 
flightpaths.     The  first  of  these  flightpaths  allow  the  PENAIR  program  to 
select  one  of  18  different  evasive  maneuvers  for  the  aircraft  upon  receiv- 
ing fire.     The  18  evasive  maneuvers  consist  of  six  basic  maneuvers  with 
three  possible  speed  changes  for  each.     The  six  basic  maneuvers  are: 
(1)  continue  flying  straight  ahead  at  the  same  altitude,  (2)  commence  a 
climb  straight  ahead  at  the  aircraft's  maximum  rate  of  ascent,  (3)  turn 
away  from  the  weapon  remaining  at  the  same  altitude,  (4)  turn  away  and 
commence  a  maximum  rate  climb,   (5)  dive  for  the  deck  and  commence  nap 
of  the  earth  (contour)  flight  straight  ahead,  or  (6)  turn  away  from  the 
weapon,  dive  for  the  deck  and  commence  nap  of  the  earth  flight.     For  each 
of  the  above  maneuvers  the  aircraft  may  decrease  speed,  remain  at  the 
same  speed  or  increase  speed.    All  of  the  above  maneuvers  are  gener- 
ated using  the  input  characteristics  of  the  aircraft  being  simulated.     It 
should  be  realized  by  the  PENAIR  user  however  that  only  one  basic  maneuver 


may  be  employed  in  one  computer  run.     The  maneuver  employed  is  deter- 
mined by  the  numerical  value  of  the  input  flag  IMARK. 

The  other  two  simulated  aircraft  flightpaths  are  preplanned.     In  the 
first  of  these,  the  flightpath  is  completely  determined  in  advance  and 
read  into  the  PENAIR  program  from  data  cards.    In  the  second,  the  x  and  y 
coordinates  of  points  along  the  flightpath  are  read  in  to  the  PENAIR  pro- 
gram.    From  these  points  the  PENAIR  program  computes  a  nap  of  the  earth 
flightpath  for  the  aircraft  over  the  simulated  terrain.     In  both  of  these 
cases,  the  PENAIR  program  makes  no  deviation  from  the  pre-planned  flight - 
path  and  the  aircraft  makes  no  use  of  the  evasive  maneuvers.    A  complete 
description  of  the  inputs  needed  to  regulate  aircraft  movement  is  included 
in  Appendix  III. 
4.     Determination  of  the  Aircraft's  Survival  Probability. 

The  major  output  of  the  PENAIR  program  is  the  probability  of  aircraft 
survival .     This  survival  probability  is  computed  as  a  function  of  weapon 
type,  range  and  time  in  view.     The  necessary  inputs  are  weapon  character- 
istics and  the  single  time-interval  kill  probabilities  for  each  weapon. 
Given  these  kill  probabilities,  the  probability  of  an  aircraft  surviving  any 
single  weapon  is  computed  as: 

P  (Survival)  =  (1  -  P. )  ( 1  -  P  J    .  .  .  (1  -  P  ) 

i  z  n 

where  P.  is  the  probability  of  kill  for  one  time  interval  at  range  i     for  a 
given  weapon.     The  total  probability  of  aircraft  survival  is  then  computed 
as  the  product  of  the  individual  probabilities  of  survival  for  each  weapon. 
(See  Flow  Chart  5.  1,  Appendix  I,  page  48  )  .     It  is  assumed  that  the  proba- 
bility of  a  kill  by  one  weapon  is  independent  of  all  other  weapons .    The 
number  of  time  intervals  the  aircraft  receives  fire  depends  on  the  time  in 
view,  minimum  and  maximum  weapon  range,  and  the  weapon  acquisition, 
firing ,  and  reloading  times  . 


Time  in  view  is  the  one  parameter  to  be  determined.  The  method  used 
by  this  simulation  to  determine  time  in  view  will  be  explained  using  Figure 
2  and  Flow  Chart  7.0  of  Appendix  I,  page  49  . 

The  aircraft  is  considered  in  view  of  a  weapon  if  every  point  on  the 
line  of  sight  WA  between  the  weapon  position  W,  and  the  aircraft  A,  is 
higher  than  the  corresponding  point  of  the  terrain.     To  determine  this  the 
x,y  coordinates  of  the  point  on  the  line  TR,  10  meters  from  the  weapon  are 
used  to  calculate  the  elevation  of  the  terrain  at  that  point.     This  elevation 
is  then  compared  with  the  height  of  the  line  of  sight  at  those  coordinates. 
If  the  terrain  is  lower,  another  point  50  meters  farther  along  the  line  TR 
is  chosen  and  the  same  comparison  is  made.    If  the  terrain  elevation  is 
below  the  line  of  sight  elevation  at  this  point  this  procedure  is  repeated 
every  50  meters.    This  process  is  continued  until  it  is  determined  that  the 
aircraft  is  in  view  of  the  weapon  or  that  the  terrain  elevation  exceeds  the 
line  of  sight.    If  the  latter  occurs  the  aircraft  is  considered  not  in  view 
of  this  weapon  at  this  time.    This  process  is  repeated  for  all  weapons  and 
then  the  aircraft  is  advanced  one  time  interval.    The  number  of  consecutive 
time  intervals  the  aircraft  remains  in  view  of  any  one  weapon  is  then  used 
to  compute  the  current  total  amount  of  time  in  view  per  weapon . 


FIGURE    2 

DETERMINATION  OF  THE  AIRCRAFT'S  TIME 
IN  VIEW  OF  A  WEAPON 


Two  other  methods  of  determining  if  the  aircraft  was  in  view  of  the 
weapon  were  tested.     Elevation  comparisons  between  line  of  sight  and 
terrain  were  made  at  2  5  meter  increments  on  the  line  of  sight  WA,  and  at 
25  meter  increments  along  TR.     No  significant  differences  were  noted  in 
the  times  of  view  computed  by  the  3  methods  mentioned.    However,  the 
method  outlined  in  the  previous  paragraph  greatly  reduces  computer  run- 
ning time  for  the  PENAIR  program. 

After  the  aircraft  has  been  in  view  of  a  weapon  for  the  required  acqui- 
sition time,  the  weapon  minimum  and  maximum  ranges  are  checked.     If 
the  aircraft  is  within  firing  range,  the  weapon  commences  firing  and  the 
probability  of  survival  is  determined  for  this  weapon  in  this  time  interval. 
The  probabilities  of  aircraft  survival  for  each  weapon  are  then  computed 
for  all  following  time  intervals  until  the  aircraft  either  goes  out  of  view 
and/or  range  or  until  the  weapon  has  expended  its  ammunition.     In  both 
cases  no  further  probabilities  are  computed  until  the  weapon  is  reloaded 
or  until  the  aircraft  comes  back  in  view  and  within  range.    It  is  assumed 
that  after  reloading  reacquisition  requires  one-half  the  acquisition  time. 
If  the  aircraft  goes  out  of  view  the  weapon  must  wait  its  full  acquisition 
time  after  the  aircraft  is  again  in  view . 
5 .     Terrain  Simulation . 

In  any  computer  model  of  a  ground  combat  environment  it  is  essential 
to  have  a  realistic  simulation  of  the  terrain  in  order  to  determine  the 
mask  angles  of  ground  weapon  positions .    To  do  this  the  simulation  must 
provide  the  elevation  of  any  point  in  the  combat  area.    Present  computer 
methods  of  terrain  simulation  either  use  very  specialized  equipment  or 
involve  prohibitive  amounts  of  computer  storage.    The  common  method  is 
to  divide  the  terrain  of  interest  into  equal  areas  and  assign  an  appropriate 
elevation  to  each.     The  assigned  elevation  might  be  the  average  elevation 
of  the  area.     In  this  method  it  is  assumed  that  all  points  in  a  given  area, 
for  the  purpose  of  simulation,  have  the  same  elevation.     This  is  a 


crude  approximation  of  terrain  features  for  large  areas .    The  approximation 
approaches  reality  as  the  size  of  the  areas  decreases.    However,  since 
the  number  of  the  areas  considered  is  inversely  proportional  to  the  size  of 
the  areas,  as  the  approximation  of   the  terrain  approaches  reality  the  com 
puter  storage  required  becomes  very  large. 

For  the  model  developed  here  the  terrain  features  are  expressed  in 
polynomial  form,  z=f(x,y)  .     The  general  problem  is  that  of  fitting  points 
in  three  dimensions  to  a  surface  represented  by  the  polynomial,  z  =  f(x,y)  . 
A  method  of  forming  this  polynomial  using  the  principle  of  least  squares, 
given  the  set  of  data  points  f  (x,y,z)  }  ,  had  been  programmed  in  FORTRAN-60 
for  the  CDC  1604  [2~\.    This  program  was  converted  to  FORTRAN-63  and  a  suc- 
cessful attempt  was  made  in  fitting  a  polynomial  to  a  fairly  rough  surface 
(mountainous  terrain.)    The  polynomials  that  best  fit  mountainous  terrain 
whose  elevations  ranged  from  800  to  1300  feet  have  approximately  400  terms 
with  the  maximum  degree  of  both  x  and  y  near  20 . 

The  program  designed  to  fit  the  polynomial  to  the  terrain  required  the 
input  points  be  from  an  equally  spaced  grid  in  the  x,y  plane.    For  this  thesis, 
areas  of  1000  meters  on  a  side  were  arbitrarily  chosen  from  maps  of  actual 
terrain.    In  order  that  the  fit  be  acceptable  on  the  edges  of  the  grid  square, 
an  area  of  at  least  1200  x  1200  meters  with  the  desired  1000  meter  square 
centered  inside  had  to  be  considered.    Consecutive  x  and  y  coordinates  of 
each  individual  contour  line  were  recorded  and  a  program  was  written  to 
produce  equally  spaced  elevations  from  the  consecutive  data  points .    How- 
ever, because  of  the  great  amount  of  storage  required,  these  two  programs 
could  only  be  incorporated  into  one  by  the  use  of  overlays.    An  overlay  is 
a  feature  of  FORTRAN-63  designed  to  increase  the  use  of  the  core  storage 
of  the  computer  by  using  magnetic  tape  as  auxiliary  storage. 

The  program  to  compute  z=  f(x,y)  from  the  consecutive  (x,y)  points 
of  actual  contour  lines  consists  of  a  main  program  and  three  overlays . 


The  sole  purpose  of  the  main  program  is  to  call  the  individual  overlays 
into  operation.     The  first  overlay,  MAPPING,  computes  the  equally 
spaced  elevation  points  using  the  (x,y)  points  from  the  contour  lines  as 
input..    The  second  overlay,   SURFIT,  by  using  the  points  computed  in  the 
first  overlay,  computes  the  polynomial  z=  f(x,y)  and  records  the  co- 
efficients on  punched  cards.     The  second  overlay  also  computes  the  ele- 
vations of  a  grid  for  use  in  the  third  overlay.     The  third  overlay,,  PLOTTER, 
generates  data  and  stores  it  on  a  magnetic  tape  suitable  for  a  CDC-160 
program  used  to  drive  a  CalComp  165  graph  plotter.     The  contour  map  drawn 
can  then  easily  be  compared  to  the  original  map  to  estimate  the  goodness 
of  fit  of  the  polynomial.     The  drawing  of  this  map  is  optional  and  is  con- 
trolled by  an  input  parameter  to  the  main  program.    A  detailed  description 
of  each  overlay  follows. 

Overlay  1  determines  the  elevation  of  the  terrain  at  equally  spaced 
points  of  a  grid  in  the  x,y  plane.     The  program  was  written  by  T.  K.  Rode- 
burg  of  the  Norwegian  Defense  Institute  for  Stanford  Research  Institute, 
Fort  Ord,  California.    This  overlay  uses  as  input  the  lists  of  the  x  and  y 
coordinates  at  successive  points  on  each  contour  line  read  from  an  en- 
larged map  of  the  actual  terrain.     The  following  steps  are  taken  by  this  over- 
lay in  the  process  of  determining  the  elevation  of  the  grid  points: 

(1)  The  point  coordinates  are  corrected  for  distortion  caused  by  un- 
equal magnification  of  the  original  map. 

(2)  Equally  spaced  vertical  lines  are  determined  and  the  coordinates 
of  the  intersection  of  the  contour  lines  and  these  vertical  lines  are  calcu- 
lated. 

(3)  The  elevation  at  these  intersections  are  calculated  by  linear 
interpolation . 

(4)  By  further  linear  interpolation,  elevations  at  equally  spaced 
points  on  the  vertical  lines  are  found,  thus  obtaining  the  elevations  of 
points  on  a  grid  in  the  x,y  plane. 
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(5)  Using  the  original  data  points,  the  above  interpolation  procedure 
is  repeated  at  equally  spaced  horizontal  lines. 

(6)  The  average  of  the  two  values  thus  obtained  is  used  as  the  ele- 
vation at  each  point  of  the  grid. 

Overlay  2  fits,  in  the  least  squares  sense,  a  polynomial  function  of 
two  independent  variables,  x  and  y,  to  values  of  a  dependent  variable,  z, 
specified  at  points  on  a  rectangular  grid  in  the  plane  of  the  independent 
variables,     That  is,  given  a  set  of  ground  elevations,  z  at  grid  coordinates 

N      M  Jim 

x  and  ye  this  overlay  finds  the  polynomial,      £      £     A      X  Y    ,  which  best 

n=o  m=o    nm 

fits  the  data.    The  use  of  orthogonal  polynomials  leads  to  a  particularly 
simple  system  of  linear  equations  rather  than  the  ill-conditioned  system 
which  arises  from  the  usual  least  squares  normal  equations.    It  also  provides 
a  measure  of  the  improvements  resulting  from  each  new  term  included  which 
further  leads,  in  this  algorithm,  to  an  automatic  selection  of  a  "best"  degree 
polynomial  function  as  determined  by  Gauss'  criterion  [1]  .    It  therefore  de- 
termines the  optimum  values  of  N  and  M  as  well  as  the  coefficients. 

This  algorithm  was  programmed  in  FORTRAN-60  for  the  CDC  1604  by 
Mrs.  Jean  Bow  of  the  USNPGS  from  an  ALGOL  program  [2]  .    Minor  modifi- 
cations were  made  to  convert  it  to  a  FORTRAN-63  overlay.    The  FORTRAN-60 
version  is  a  very  general  program  whereas  the  overlay  has  had  many  of  the 
generalities  removed  for  ease  of  handling  inputs  which  in  turn  reduced  the 
computer  running  time  significantly.    A  FORTRAN  listing  of  this  program 
appears  in  Appendix  IV, 

The  following  steps  are  taken  by  Overlay  2  in  the  process  of  determin- 
ing an  approximating  polynomial: 

(1)  the  variables  are  normalized  by  subtracting  off  their  means, 

(2)  the  orthogonal  polynomials  are  evaluated  , 

(3)  the  contribution  of  each  orthogonal  polynomial  to  the  minimization 
of  the  residuals  is  evaluated, 


(4)  Gauss'  criterion  is  applied  to  determine  the  degree  of  poly- 
nomial which  yields  the  closest  fit  to  the  given  data. 

(5)  the  orthogonal  polynomial  coefficients  are  evaluated, 

(6)  the  dependent  variables  are  evaluated  using  the  approximating 
polynomial,  and 

(7)  the  terrain  elevations  utilizing  a  50  meter  grid  are  evaluated  for 
input  into  Overlay  3 . 

Overlay  3  generates  lists  of  x  and  y  coordinates  of  successive 
points  on  each  contour  line  and  stores  these  coordinates  on  an  output 
tape  in  a  form  suitable  for  a  CDC-160  program  used  to  drive  a  CalComp  165 
plotter.  This  tape  is  used  as  an  input  tape  for  the  plotter  which  then  recon- 
structs the  terrain  by  connecting  the  points  with  straight  line  segments  . 

The  FORTRAN-63  program  used  in  Overlay  3  was  originally  developed 
for  use  with  the  IBM  7090  system  and  an  x,y  plotter,  to  give  a  complete 
traditionally  styled  contour  map.     This  program  was  converted  to  FORTRAN- 
63  by  Mr.  R.  R.  Hillery  of  the  USNPGS,  for  use  on  the  CDC  1604.  Subrou- 
tines used  to  compute  the  points  to  be  plotted  were  provided  by  Dr,   M.  O. 
Dayhoff  of  the  National  Biomedical  Research  Foundation,  Silver  Spring, 
Maryland  [3]  .  Subroutines  used  internally  to  write  the  output  tape  were 
programmed  by  Mr.   J.  R.  Ward  of  the  USNPGS  for  use  with  a  CalComp 
165  graph-plotting  routine,  DRAW. 

Figure  3  is  an  actual  contour  map  from  which  approximately  1200  input 
points  have  been  read.    A  polynomial  of  400  terms  was  fitted  and  the  ele- 
vations of  a  grid  (spacing  of  50  meters)  were  computed  to  plot  Figure  4. 

The  coefficients  of  this  polynomial  ranged  from  904.18  for  the  constant 

-94  19   19 

term  to  -4.593  x  10  for  the  x    y      term. 

An  analysis  of  a  table  of  differences  showed  the  maximum  difference 

of  625  actual  and  computed  elevations  to  be  31  feet  in  an  elevation  range 

of  900  to  1200  feet. 
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FIGURE  3 

ACTUAL  CONTOUR  MAP  OF  GRID  SQUARE  7174 

U.S.  ARMY  HUNTER- LIGGETT  MILITARY 
RESERVATION 


FIGURE  4 

COMPUTER  SIMULATION  OF  GRID  SQUARE  7174 
U.  So  ARMY  HUNTER-LIGGETT  MILITARY 
RESERVATION 


11 


FIGURE  5 

REPRESENTATIVE  FIT  AT  THE  JUNCTION 
OF  TWO  GRID  SQUARES 
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Figure  5  shows  a  representative  fit  at  the  junction  of  two  of  the  grid 
squares  approximated  by  separate  polynomials.    The  accuracy  with  which 
the  grid  squares  fit  together  is  determined  by  the  accuracy  of  the  input 
data. 
6.    Analysis  of  the  Model 

PENAIR  has  been  checked  out  with  the  use  of  a  simple  scenario.    A 
3000  x  4000  meter  mountainous  section  of  the  U.  S.  Army  Military  Reser- 
vation was  selected  for  simulation.     Three  10 -men  enemy  squads  armed 
with  M-14  rifles  were  situated  as  illustrated  in  Figure  6.    A  UH-1B  heli- 
copter was  assigned  the  mission  of  delivering  supplies  from  an  airfield 
at  position  a  to  an  outpost  at  position  b    ,     Two  test  runs    were  made  with 
altitudes  of  the  helicopter  of  1000  feet  and  500  feet  above  the  terrain.    In 
both  tests  all  18  evasive  maneuvers  were  checked.     Some  representative 
flightpaths  generated  by  the  program  are  illustrated  in  Figures  6 ,  7 ,  8 ,  and 
9 .    The  survival  probabilities  of  the  aircraft  on  the  various  attempts  were 
computed  and  are  shown  in  Tables  1  and  2 .    All  of  the  computed  proba- 
bilities are  realistic  and  are,  in  a  comparative  sense,  what  one  would 
expect. 

A  sensitivity  analysis  was  then  conducted  on  the  rate  of  climb  and  des- 
cent and  changes  in  aircraft  speeds .    A  helicopter  was  simulated  flying  at 
100  knots  with  its  maximum  rate  of  climb  and  descent  taking  the  values   of 
500,1000  and  2000  feet  per  minute.    The  effects  of  these  three  rates  on 
nap  of  earth  flight  are  illustrated  in  Figure  10  .    It  should  be  noted  that 
Figure  10  illustrates  how  the  program  generates  the  flightpath  to  overshoot 
the  crest  of  a  hill  before  descending  just  as  a  pilot  tends  to  do.    Figure  10 
also  shows,  as  should  be  expected,  the  greater  the  rate  of  climb  and  des- 
cent, the  closer  the  aircraft  remains  to  the  earth.    Figure  11  illustrates  an 
aircraft  flown  with  a  maximum  rate  of  climb  and  descent  of  1500  feet  per 
minute  and  speeds  of  100,  200  and  300  knots.     The  resulting  flightpaths  also 
appear  realistic. 
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The  flightpath  computed  by  PENAIR  illustrated  in  Figure  12  is  of  a 
helicopter  making  a  level  turn  at  70,  90  and  110  knots,     Here  it  is  noted, 
for  the  computed  flightpath,  the  radius  of  turn  increases  as  the  airspeed 
increases . 

The  last  analysis  made  was  the  comparison  of  computed  results  with 
actual  results.     Figure  13  is  the  simulated  flightpath  and  terrain  of  an 
actual  flightpath  flown  by  a  UH-1B  over  an  actual  terrain      The  44  weapon 
positions  illustrated  were  actually  manned  by  U ,  S,   Army  personnel.  The 
observed  time  the  helicopter  was  in  view  of  each  of  these  positions  is 
listed  in  Table  3.     The  times  in  view  for  the  simulated  terrain  and  flight- 
path  were  computed  and  are  also  listed  in  Table  3.    It  is  believed  that  with 
better  terrain  data  an  even  closer  comparison  could  be  obtained. 
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FIGURE  6 


A  STRAIGHT  AND  LEVEL  FLIGHTPATH 
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Terrain  Profile 
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FIGURE  7 


A  FLIGHTPATH  OF  AN  AIRCRAFT  CLIMBING  STRAIGHT  AHEAD 
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FIGURE  8 


A  FLIGHTPATH  OF  AN  AIRCRAFT  TURNING  AWAY  FROM  A  FIRING  WEAPON 
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FIGURE  9 


A  NAP  OF  THE  EARTH  FLIGHTPATH 
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Summary  of  Survival  Probability  of  A  UH-1B  Helicopter  Fired  on  by  Three 
10-Men  Squads  of  M-14  Rifles 


Initial  Speed  90  Kts 


Initial  Altitude  650  Meters 


Evasive  Maneuver  Evasive 

(Commenced  after  receiving  fire)         Speed 


Survival  Probabilities 
Wpn  A     Wpn  B      Wpn  C      Total 


Fly  straight  and  level 


70  .8576        .8078        1.0000     .6928 

90  .8723         .8649        1.0000     .7545 

110  .8723        .8723        1.0000      .7610 


Climb  straight  ahead 


70  .8723        .8723        1.0000     .7610 

90  .8723        .8723        1.0000     .7610 

110  .8723        .8723        1.0000     .7610 


Turn  away,  level 


70  .9026      1.0000        1.0000     -9026 

90  .9026      1.0000        1.0000      -9026 

110  .9026      1.0000        1.0000      .9026 


Turn  away  and  climb 


70  .9026      1.0000 

90  .9026      1.0000 

.110  .9026      1.0000 


1.0000  .9026 
1.0000  .9026 
1.0000     .9026 


Dive  and  fly  NOE 


70  .8359        .7252        1.0000      .6062 

90  .8723        .8031        1.0000      .7005 

110  .8723        .8242        1.0000     .7190 


Turn  away  and  fly  NOE 


70  .8798      1.0000        1.0000     .8798 

90  .8874      1.0000        1.0000     .8874 

110  .8874      1.0000        1.0000      .8874 


Table  1 
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Summary  of  Survival  Probabilities  of  UH-1B  Helicopter  Fired  on  by  Three 

10-Men  Squads  Armed  with  M-14  Rifles 


Initial  Speed  90  Kts 


Initial  Altitude  450  Meters 


Evasive  Maneuver 
(Commenced  after  receiving  fire) 


Evasive  Survival  Probabilities 

Speed    Wpn  A      Wpn  B       Wpn  C     Total 


Fly  straight  and  level 


70 

90 
110 


Values  not  determined 
at  these  speeds 

8309  .7788     1.0000 


6471 


Climb  straight  ahead 


70 

90 
110 


Values  not  determined 
at  these  speeds 

8377        .8242        1.0000       .6905 


Turn  away,  level 


70         .8874      1.0000        1.0000       .8874 

90  .8874      1.0000        1.0000       .8874 

110         .8874      1.0000        1.0000       .8874 


Turn  away  and  climb 


70         .8798      1.0000        1.0000       .8798 

90         .8798      1.0000        1.0000       .8798 

110         .8798      1.0000        1.0000       .8798 


Dive  and  fly  NOE 


70         .7558         .9501         1.0000       .7181 

90         .8027         .9664        1.0000       .7758 

110         .3306      1.0000        1.0000       .8306 


Turn  away  and  fly  NOE 


70         .8798      1.0000        1.0000       .8798 

90         .8798      1.0000        1.0000       .8798 

110         .8798      1.0000        1.0000       .8798 


Table  2 
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FIGURE  10 

COMPARISON  OF  RATES  OF  CLIMBS  AND  DESCENTS  IN  A  NAP 
OF  THE  EARTH  FLIGHT  AT  A  CONSTANT  SPEED  OF  100  KTS. 
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COMPARISON  OF  THE  AIRCRAFT  S  RADIUS    OF  TURN 
AT  3  AIRSPEEDS,  70,  90  and  110  KTS. 
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COMPUTER  SIMULATION  OF  AN  ACTUAL  FLIGHTPATH 
FLOWN  BY  A  UH- -IB  HELICOPTER 
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Comparison  of  Actual  and  Computed  Time  Intervals  in  View  of  44  Weapon  Sites 

Time  Interval  =    .6  sec. 


Wpn. 

Time 

in    View 

Wpn 

Time    ] 

n    View 

Wpn 

Time 

In    View 

No 

Actual 

Computed 

No 

Actual 

Computed 

No 

Actua] 

Computed 

1 

0 

0 

16 

50 

56 

31 

0 

0 

2 

0 

9 

17 

0 

0 

32 

22 

39 

3 

30 

30 

18 

43 

47 

33 

0 

0 

4 

0 

0 

19 

63 

108 

34 

89 

88 

5 

0 

0 

20 

4 

30 

35 

93 

91 

6 

0 

59 

21 

17 

74 

36 

69 

76 

7 

0 

0 

22 

0 

72 

37 

66 

58 

8 

0 

40 

23 

0 

0 

38 

49 

64 

9 

0 

0 

24 

0 

0 

39 

16 

36 

10 

0 

17 

25 

59 

63 

40 

0 

16 

11 

9 

2 

26 

0 

34 

41 

8 

8 

12 

0 

0 

27 

63 

64 

42 

56 

56 

13 

26 

39 

28 

43 

45 

43 

19 

15 

14 

40 

63 

29 

0 

9 

44 

103 

102 

15 

55 

73 

30 

0 

0 

Table    3 
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7.     Conclusions  and  recommendations. 

The  simulation  of  the  terrain  obtained  by  this  program  has  proven  to  be 
extremely  realistic.  (See  Figure  4,  Section  5,  page  11   )  .     The  junctions  of 
squares  approximated  by  separate  polynomials  are  as  good  as  the  data  used 
for  the  approximation.   (See  Figure  5,  Section  5,  page  12  )  .    This  program  was 
written  in  its  present  form  because  of  the  availability  of  input  data  to  the 
author  in  the  format  described  in  Appendix  V.     If  the  actual  elevations  could 
be  obtained  using  an  equally  spaced  grid  for  the  x  and  y    coordinates,  the 
first  overlay  could  be  eliminated      It  is  felt  that  this  is  feasible  and  would 
make  the  preparation  of  the  input  cards  much  simpler. 

Among  the  uses  of  the  simulation  of  the  terrain  by  this  method  include 
ground  to  ground  combat  and  ground  to  air  combat  environments .    Mask  angles , 
terrain  profiles  and  actual  distances  between  points  in  a  mountainous  terrain 
are  all  easily  determined.     The  uses  of  the  PENAIR  program  include  the  deter- 
mination of  optimal  flightpaths  across  hostile  terrain,  the  survival  proba- 
bility and  optimal  maneuvers  to  perform  if  the  aircraft  is  fired  upon. 

The  results  in  Section  6  all  tend  to  illustrate  the  realism  of  the  aircraft 
penetration  model.    Whether  the  actual  kill  probabilities  used  are  exact  is 
not  extremely  important  since  an  appropriate  defensive  maneuver  for  an  un- 
armed aircraft  to  initiate  upon  receiving  ground  fire  may  be  obtained  by  ob- 
serving the  relative  survival  probabilities  of  several  maneuvers. 

It  is  believed  that  with  additional  programming  this  model  could  be  ex- 
tended to  include  formations  of  two  or  more  aircraft ,  thus  greatly  increasing 
its  usefulness.    In  addition,  since  the  terrain  is  simulated  by  mathematical 
functions  and  since  the  line  of  sight  between  aircraft  and  weapon  may  be 
described  by  a  linear  relationship,  it  appears  possible  to  solve  the  two  rela- 
tions for  a  point  of  intersection.    Upon  accomplishing     this  programming,  it 
is  believed  the  computer  running  time  would  be  reduced  significantly. 
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APPENDIX  I 
Flow  Charts  of  the  PENAIR  Program 

In  the  flow  charts  of  this  appendix,  the  following  conventions 
are  used: 


This  symbol  is  used  for  a  READ 
statement. 


This  symbol  is  used  for  a  PRINT 
statement. 


A  rectangle  indicates  one  or  more 
FORTRAN  statements  .    A  number  such 
as  4.1.8  indicates  the  statement 
in  this  rectangle  is  further  explained 
by  Flow  Chart  4.1.8  in  this  appendix. 
A  number  such  as  420  is  the  statement 
number  as  it  appears  in  the  FORTRAN 
listing  in  Appendix  II. 


This  symbol  is  used  for  an  IF  statement, 


© 


This  symbol  indicates  the  logic  of  this 
Flow  Chart  is  continued  at  the  beginning 
of  Flow  Chart  5.1  of  this  appendix. 
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A 


This  symbol  indicates  a  connector  between 
two  parts  of  the  same  Flow  Chart  of  this 
appendix. 
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This  symbol  indicates  the  logic  of  this 
Flow  Chart  is  continued  where  the  same 
symbol  appears  in  Flow  Chart  4.1  of  this 
appendix. 

This  is  the  symbol  for  a  CONTINUE  state- 
ment whose  statement  number  is  10. 

This  is  the  symbol  for  a  computed 
GO  TO  statement. 


This  symbol  is  used  for  the  statement 


DO    10    1=  1,  12 


BE&IN 


This  symbol  indicates  the  beginning  or 
entry  point  for  the  Flow  Chart . 


END 


This  symbol  indicates  the  end  or  exit 
point  for  the  Flow  Chart. 
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APPENDIX  II 

FORTRAN  Listing  of  the  PENAIR  Program 

PENAIR  consists  of  the  main  program  and  four  subroutines.     The 
variable  names  used  in  the  PENAIR  program  are  defined  in  the  comments 
appearing  in  the  FORTRAN  listing.    The  variable  names  of  input  and  out- 
put parameters  are  explained  in  Appendix  III. 
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Appendix  III 
Inputs  and  Outputs  of  PENAIR 
Data  cards  for  the  PENAIR  program  should  be  prepared  in  the  following  manner 


Data 

Number 

Variable 

Columns 

Description 

Group 

of  Cards 

Name 

of    Card 

of    Variable 

1 

1 

NRUN 

1-10 

A  right  justified  integer, 
the  number  of  runs  to  be 
performed . 

KMAX 


1-10 


XMIN 


YMIN 


XMAX 


11-20 


21-30 


31-40 


A  right  justified  integer, 
the  number  of  squares  used 
to  approximate  the  terrain, 
KMAX  £  12 « 

A  right  justified  floating 
point  number,  the  x  co- 
ordinate of  the  left  edge  of 
the  terrain .    All  coordinates 
are  in  meters ,  that  is  the 
distance  between  the  two 
x  coordinates  72455  and 
72454  is  one  meter. 
A  right  justified  floating 
point  number,  the  y  co- 
ordinate of  the  bottom  edge 
of  the  terrain . 

The  x  coordinate  of  the  right 
edge  of  the  terrain. 
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Date 

Group 


Number 
of  Cards 


Variable 
Name 

YMAX 


TI 


TPD 


X2 
Y2 


GMAX 


RC 


Columns 
of  Card 

41-50 


51-60 


61-70 


X 

1-10 

Y 

11-20 

Z 

21-30 

31-40 
41-50 


1-10 


11-20 


Description 
of  Variable 

The  y  coordinate  of  the 
top  edge  of  the  terrain. 
A  right  justified  floating 
point  number,  the  time 
interval  of  the  game 
(seconds)  ,  TI  >  0  . 
A  right  justified  floating 
point  number,  the  dimen- 
sion of  the  map  square 
approximated  by  one  poly- 
nomial (meters)  . 

Right  justified  floating 
point  numbers ,  the  begin- 
ning coordinates  of  the 
aircraft.     If  IMARK  =  3,4,  or 
6,.  y  must  equal  YMIN. 
The  x  and  y  coordinates  of 
the  desired  destination  of 
the    aircraft .     If  IMARK  =  7 
or  8  „  this  card  may  be  blank . 

A  right  justified  floating 
point  number,  the  maximum 
positive  G  force  to  be  used 
by  the  aircraft. 
A  right  justified  floating 
point  number,  the  aircraft's 
maximum  rate  of  climb 
(feet/minute) 


77 


Data 
Group 


Number 
of  Cards 


Variable 
Name 

RD 


GS 


Columns 
of  Card 

21-30 


31-40 


Description 
of  Variable 

A  right  justified  floating 
point  number,  the  aircraft  s 
maximum  rate  of  descent 
(feet/minute)  . 
A  right  justified  floating 
point  number,  the  original 
airspeed  of  the  aircraft. 


IMAX 


1-10  A  right  justified  integer, 

the  number  of  types  of 
weapons ,  IMAX  £  10  .   . 


WPNAWPNB     1-16 


JMAX 


TA 


FT 


21-25 


26-30 


31-35 


The  name  of  the  weapon  „ 
e.g.  ,  50  Cal.  MG 
A  right  justified  integer, 
the  number  of  weapons  of 
this  type,    JMAX  <  20. 
A  right  justified  integer, 
the  time  required  to  acquire 
the  target  by  this  weapon 
(seconds)  ,  TA  £  0  . 
A  right  justified  integer, 
the  firing  time  of  this 
weapon  before  reloading 
is  required  (seconds)  , 
FT>  0. 
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Data 
Group 


Number 
of  Cards 


Variable 
Name 

RLT 


RMAX 


RMIN 


Columns 
of  Card 

36-40 


41-50 


51-60 


Description 
of  Variable 

A  right  justified  integer, 
the  time  to  reload  this 
weapon  (seconds)  RLT  £  0  . 
A  right  justified  floating 
point  number,  the  maximum 
range  of  this  weapon(meters) 
RMAX  >  0  . 

A  right  justified  floating 
point  number,  the  minimum 
range  of  this  weapon(meters) 
RMIN  ;>  0  . 


1-8  The  probability  of  kill  by 

9-16  this  weapon  in  one  time 

interval  at  ranges  of 
72-80  100,  200,  300,    ...,  1000 

meters . 
1-8  The  second  card  of  this 

9-16  group  contains  the  kill 

probabilities  at  ranges  of 
72-80  1100  to  2000  meters  in  100 

meter  intervals . 


JMAX 


NW 


1-2 


xw 

11-20 

YW 

21-30 

ZW 

31-40 

The  number  assigned  to  the 
weapon,    1  ^  NW  £  99  . 

The  x,  y  and  z  coordinates 
of  this  weapon,  ZW  may  be 
blank  since  PENAIR  will  com- 
pute the  weapon's  elevation. 
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NOTE:    Data  groups  6 ,  7  and  8  should  be  repeated  IMAX  times . 


Data 
Group 


Number 
of  Cards 


Variable 
Name 

MMAX 


NMAX 


Columns 
of  Card 

1-10 


11-20 


Description 
of  Variable 

A  right  justified  integer, 
the  degree  of  x  in  the  poly- 
nomial that  approximates 
the  first  map  square. 
A  right  justified  integer, 
the  degree  of  y  in  the 
polynomial  that  approximates 
the  first  map  square 


10 


As  many  cards      PHI 
as  necessary 


1-80  The  coefficients  of  the 

first  polynomial  appear  on 
this  and  the  following  cards 
of  this  data  group  in  a 
4E2010  field  as  obtained 
from  the  terrain  simulation 
program . 


NOTE:  Data  groups  9  and  10  should  be  repeated  KMAX  times . 


11 


AFIELD  1-16  These  four  cards  contain 

the  four  lines  of  the  title 
to  be  printed  at  the  top  of 
the  output. 
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Data  Number  Variable      Columns  Description 

Group  of  Cards  Name       of  Cards  of  Variable 

12  1  IMARK  10  The-value  of  this  flag,  deter- 

mines the  flightpath  to  be  flown. 
IMARK  =  1,  upon  receiving  fire, 
fly  straight  and  level , 
=  2,  climb  straight  ahead, 
=  3,  turn  away,  level, 
=  Ac  turn  away  and  climb, 
=  5  0  dive  and  fly  nap  of 

the  earth, 
=  6 ,  turn  away  and  fly 
nap  of  the  earth, 
=  7 ,  the  flight  path  is 
completely  pre- 
determined and  read 
in  as  data  group  13, 
=  8 ,  the  x  and  y  coordi- 
nates are  read  in 
as  data  group  13, 
and  a  nap  of  the  earth 
flightpath  is  computed, 
IPUNCH         20  IPUNCH  =  1 ,   the  flightpath  will 

will  be  punched  on 

cards . 
=    0,  the  flightpath  will 

not  be  punched. 


Data 

Number 

Variable 

Column 

Description 

Group 

of  Cards 

Name 

of  Card 

of  Variable 

GSNEW 

21-30 

A  i 

ight  justified 

point  number,  the  ground 
speed  to  be  attained  by  the 
aircraft  after  receiving  fire . 


NOTE:    Data  groups  11  and  12  will  be  repeated  NRUN  times. 


13 


As  many 
cards  as 

NPATH 

6-7 

necessary 
but  not 
greater 
than  300. 

XT 

YT 

18-26 
27-35 

ZH 

36-44 

Any  positive  integer,  used 
to  identify  the  flightpath. 

Right  justified  floating 
point  numbers ,  each  of  the 
cards  of  this  group  contain 
x,  y,  and  z  coordinates  of 
the  aircraft  inconsecutive 
time  intervals . 


NOTE:    A  blank  card  should  be  inserted  at  the  end  of  this  data  group. 
This  data  group  will  be  used  only  if  IMARK  =  7  or  8 .    If 
IMARK  =  8,  ZH  may  be  blank.     Data  group  13  should  follow 
its  corresponding  data  groups  11  and  12. 
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Output  of  the  PENAIR  Program 

A  sample  output  from  the  PENAIR  program  is  illustrated  in  Figures 
14  and  15      The  input  parameters  of  the  weapons  are  printed  with  the 
weapon  positions  and  elevations.     The  second  page  of  the  output  (Figure 
15)  is  a  summary  of  time  in  view  and  aircraft  survival  probability.  The 
time  in  view  listed  here  is  in  number  of  time  intervals  in  view.     If  the 
time  interval  used  in  this  example  had  been  0.5  second,  the  aircraft 
was  in  view  of  weapon  number  49,  located  at  (70080,  76940,  339)  ,  for 
28  seconds  (56x0.5).     The  survival  probability  of  the  aircraft  from  this 
weapon  was   .9991. 
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APPENDIX  IV 

FORTRAN  Listing  of  the  TERRAIN  Program 

TERRAIN  consists  of  a  main  program  and  three  overlays.    The  control 
cards  for  the  COOP-Monitor  of  the  CDC  1604  are  included  in  this  listing 
and  are  explained  in  Appendix  V.    The  variable  names  used  in  the  TERRAIN 
program  are  defined  in  the  comments  appearing  in  the  FORTRAN  listing. 
The  variable  names  of  input  and  output  parameters  are  explained  in  Appendix 
V. 
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APPENDIX  V 

Inputs  and  Outputs  of  the  Terrain  Simulation  Program 

The  following  information  is  presented  for  potential  users  of  this  pro- 
gram: 

A.  Tape  units 

(1)  The  COOP  control  cards  are  shown  in  the  FORTRAN  listing  (See 
Appendix  IV.)    A  dash  (-)  in  column  one  indicates  a  7,9  punch,  while  an 
open  parenthesis  (()  indicates  an  11,0,7,9  punch , 

(2)  Logical  tape  unit  49  is  the  graph  tape  and  45  refers  to  the  type- 
writer . 

(3)  Scratch  tapes  56  and  57  are  necessary  and  therefore  a  total  of 
eight  tape  units  are  required 

B.  Input  Formats 

The  preparation  of  the  input  data  cards  is  explained  on  the  following  pages 
of  this  appendix. 
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Data  Number       Variable 

•  Group  of  Cards        Name 

1  2  ITITLE 


NL 


NL 


II 


LABL(ir) 


NMAX 


MMAX 


Columns 
of  Card 

1-48 


1-5 


1-5 


7-10 


1-5 


6-10 


NEWCTL  11-15 
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Description 
of  Variable 

These  two  cards  contain  the 
two  line  title  desired  for 
the  first  page  of  the  output 
and  the  bottom  of  the  plotted 
map. 

A  right  justified  integer, 
the  number  of  contour  lines 
requiring  a  label. 

A  right  justified  integer, 

the  number  of  the  contour 

level  to  be  labeled, 

1  5  II*   50 

A  left  justified  integer, 

the  elevation  of  contour 

level  II. 

A  right  justified  integer, 
the  maximum  degree  of  x  to 
be  considered 
A  right  justified  integer, 
the  maximum  degree  of  y  to 
be  considered. 
A  right  justified  integer, 
NEWCTL  =  0 ,    no  contour  map 
is  desired. 
=  1,      an  output  tape 

will  be  generated 
for  the  CalComp 
165  plotter. 


Data  Number         Variable         Columns 

Group  of  Cards  Name  of  Cards 


S 


1-10 


11-20 


Description 
of  Variable 

A  right  justified  floating 
point  number,  the  spacing 
desired  between  grid  lines  . 
S  must  be  chosen  so  that 
neither  the  number  of  vertical 
nor  horizontal  grid  lines 
exceed  45. 

A  right  justified  floating 
point  number,  the  linear 
dimension  of  the  map  square. 


XTMIN         1-10 


XTMAX       11-20 


CND  21-30 


A  right  justified  floating 
point  number,  the  elevation 
of  the  lowest  contour  in 
the  map  square. 
A  right  justified  floating 
point  number,  the  maximum 
elevation  in  the  map  square. 
A  right  justified  floating 
point  number,  the  contour 
interval .  The  first  contour 
is  the  value  C  =  XTMIN. 


GX 


1-2 


The  first  two  digits  of  a 
five  digit  floating  point 
number,  GX  *  L000  =  the  x 
coordinate  of  the  lower 
left  corner  of  the  map  square 
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Data 
Group 


Number 
of  Cards 


Variable 
Name 

GY 


XP(1) 


YP(1) 


XP(2) 
YP(2) 
XP(3) 
XP(3) 
XP(4) 
XP(4) 


Columns 
of  Cards 

3-4 


Description 
of  Variatlon 


The  first  two  digits  of  a  five 
digit  floating  point  number  , 
GY  *  1000  =    the  y  coordinate 
of  the  lower  left  comer  of  the 
map  square. 

11-13         A  right  justified  floating  point 
number,  the    x    coordinate  of 
the  lower  left  corner  of  the  map 
square  measured  on  an  arbitrary- 
scale.    XP  (1)  usually  should 
equal  0000. 

16-20       The  y  coordinate  of  the  lower  left 
comer  measured  on  an  arbitrary 
scale.    YP(1)  usually  should  equal 
0000. 

21-25       The  x  coordinate  of  the  lower 
right  corner, 

26-30       The  y  coordinate  of  the  lower 
right  corner. 

31-35       The  x  coordinate  of  the  upper 
right  corner. 

36-40       The  y  coordinate  of  the  upper 
right  comer. 

41-45       The  x  coordinate  of  the  upper 
left  comer. 

46-50      The  y  coordinate  of  the  upper 
left  comer. 
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Data 
Group 

Number 
of  Cards 

Variable 
Name 

Columns 
of  Card 

8 

As  many 
cards  as 
necessary 

1-4 

FZ 


NC 


X(l) 


Y(l) 


5-8 


9-LO 


11-15 


16-20 


21-80 


Description 
of  Variable 

Any  identifying  cede  desired 
Col.  1-4  are  not  read. 
A  right  justified  floating 
point  number,  FZ  *  10   = 
the  elevation  (feet)  whose 
coordinates  appear  on  this 
card. 

A  sequence  number  starting 
at  01  and  increasing  by  one 
on  each  card  following  whose 
coordinates  are  those  of  the 
same  contour  line. 
A  right  justified  floating 
point  number,  the  x  coordi- 
nate of  a  point  on  the  contour 
line. 

A  right  justified  floating 
point  number,  the  coordinate 
of  a  point  on  the  contour  line  „ 

F  5 . 0  fields  of  additional 
points  on  the  contour  line 
in  the  same  manner  as 
Col.  11-20. 


NOTE:  One  blank  card  is  inserted  after  this  group  of  data  cards  (eighth 
data  group)  for  each  contour  line.  One  additional  blank  card  is 
inserted  if  a  new  seventh  data  group  is  read. 
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Data 
Group 


Number 
of  Cards 

1 


Variable 
Name 


Columns 
of  Card 

7-9 


Description 
of  Variable 

The  figure  999  should  appear 
in  these  colums  to  signify- 
that  this  is  the  last  card 
of  the  data  deck.. 


NOTE:    There  should  be  no  blank  cards  between  the  last  eighth  data  group 
and  the  ninth  data  group. 
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Example: 

It  is  desired  that  grid  square  717  5  of  the  U.   S    Army  Hunter-Liggett 
Military  Reservation  be  approximated  by  a  polynomial.     Figure  17  illus- 
trates the  data  cards  necessary 
The  words 

HLMR  CLOSE  TERRAIN 
GRID  SQUARE  7175 
will  appear  at  the  top  of  the  first  page  of  the  printed  output  and  at  the 
bottom  of  the  plotted  map      The  minimum  and  maximum  elevation  appear- 
ing in  this  square  are  900  and  1200  feet  respectively.     The  contour  interval 
is  20  feet.     The  maximum  degree  of  X  and  Y  to  be  considered  is  20  .     The 
grid  interval  will  be  50  meters  and  the  square  is  1000  meters  on  a  side. 
Figure  16  illustrates  how  the  original  map  was  enlarged  to  read  the 
contour  coordinates  . 


(0012,3485) 


(0000,  0000) 


(3455,   3440) 


(3452,  0008) 


FIGURE  16 


DISTORTION  CAUSED  BY  PROJECTION 
OF  A  MAP  GRID  SQUARE 

Distortion  caused  by  projection  is  apparent  in  Figure  16  and  is  cor- 
rected by  the  first  overlay.  The  corner  coordinates  are  arbitrary  and  all 
contour  coordinates  are  relative  to  them.  If  the  points  are  measured  by 
hand,  the  actual  corner  coordinates  may  be  used  and  in  this  case  would 
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FIGURE  17 
A  SAMPLE  DATA  DECK  FOR  THE  TERRAIN  PROGRAM 
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be  (71000,   75000),   (72000,   75000-)     ,(72000,   76000),  (71000.    76000).   The 
contour  coordinates  would  then  be  the  actual  terrain  coordinates. 

In  this  example  the  grid  square  was  reprojected  a  second  time  with  new 
corner  coordinates  of  (0000,  0000),  (3444,-0020),  (3441,  3467),  (-0004   3475) 

Terrain  has  three  types  of  output ,     a)  printed ,    b)  punched  cards  ,  and 
c)  a  plotted  map  (See  Figure  18)  .     The  printed  output  consists  of    1)  the  x  and 
y  coordinates  of  the  grid  to  be  used,  2)  the  input  elevations  at  points  on  this 
grid  as  interpolated  by  Overlay  1,3)  the  elevations  of  the  same  points  as 
computed  using  the  approximating  polynomial,  4)  a  table  of  differences  be- 
tween the  actual  and  computed  elevations  of  the  grid,  5)  the  computed  ele- 
vations to  be  used  in  Overlay  3,  6)  the  coefficients  of  the  approximating 
polynomial,  and  7)  the  actual  mean  elevation  squared,  computed  mean 
elevation  squared  and  the  maximum  differences  between  the  actual  and  com- 
puted elevations.    Although  accuracy  is  not  determined  in  general,  if  the 
actual  mean  elevation  squared  and  computed  mean  elevation  squared  are 
equal,  then  computation  is  exact.    In  practice  they  will  not  be  equal  due 
to  the  imprecise  nature  of  calculation.    A  wide  discrepancy  indicates  ex- 
cessive errors  in  calculation. 

The  first  card  of  the  punched  card  output  contains  the  degree  of  x  and 
y  in  a  2110  format.     The  remaining  cards  contain  the  coefficients  of  the 
approximating  polynomial  in  a  4E20.10  format. 
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FIGURE  18 

A  SAMPLE  OF  A  PLOTTED  MAP 
FOR  THE  TERRAIN  PROGRAM 
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A  computer  simulation  of  an  aircraft  pen 
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